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We give a status report on the determination of a set of parton distributions based on 
neural networks. In particular, we summarize the determination of the nonsinglet quark 
distribution up to NNLO, we compare it with results obtained using other approaches, 
and we discuss its use for a determination of a s . 

1 Introduction 

The LHC will require an approach to the search for new physics based on the precision 
techniques which are customary at lepton machines [TJ [5] . This has recently led to significant 
progress in the determination of parton distribution functions (PDFs) of the nucleoli. The 
main recent development has been the availability of sets of PDFs with an estimate of the 
associated uncertainty [31 21 [5]. However, the standard approach to the determination of 
the uncertainty on parton distributions has several weaknesses, such as the lack of control 
on the bias due choice of a parametrization and, more in general, the difficulty in giving a 
consistent statistical interpretation to the quoted uncertainties. 

These problems have stimulated various new approaches to the determination of PDFs [7J , 
in particular the neural network approach, first proposed in Ref. [BJ. The basic idea is to 
combine a Monte Carlo sampling of the probability measure on the space of functions that 
one is trying to determine [7J with the use of neural networks as universal unbiased interpo- 
lating functions. In Rcfs. [BJ[5] this strategy was successfully applied to a somewhat simpler 
problem, namely, the construction of a parametrization of existing data on the DIS struc- 
ture function i*2(^, Q 2 ) of the proton and neutron. The method was proven to be fast and 
robust, to be amenable to detailed statistical studies, and to be in many respects superior 
to conventional paramctrizations of structure functions based on a fixed functional form. 

The determination of a parton set involves the significant complication of having to go 
from one or more physical observables to a set of parton distributions. Recently [9] most of 
the technical complications required for the construction of a neural parton set have been 
tackled and solved in the process of constructing a determination of the quark isotriplet 
parton distribution. This work will be reviewed here. Also, based on this work, we will 
present preliminary results on the determination of a s and a determination of the variation 
in x 2 which corresponds to a one-sigma variation of the underlying parton distributions. 
Work to apply the techniques of [9] to the singlet sector is at an advanced stage [TU] , 
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2 Determination of the nonsinglet quark distribution 



The first application of the neural network approach to parton distributions, a determination 
of the NS parton distribution qnsix, Qq) = (u + u — d + d)(x, Qq) from the DIS structure 
function data of the NMC and BCDMS collaborations, was presented in Ref. [9]. Results 
for this PDF were obtained at LO, NLO, NNLO for different values of a s (M§). 

In Ref. [9] we have implemented a new fast and efficient method for solving the evolution 
equations up to NNLO. This method combines the advantages of x— space and TV— space 
evolution codes: an x dependent Green function (evolution factor) is determined by inverse 
Mcllin transformation of the exact TV-space expression and stored. Evolution of PDFs is 
then performed by convoluting this Green function with any given boundary condition. The 
accuracy of this method has been benchmarked up to NNLO with the help of the tables of 
Refs. [BE]. 

Also, we have implemented a criterion to determine the convergence of the fitting proce- 
dure in a way which is free of bias related to the choice of parametrization. To this purpose, 
the dataset is randomly divided into two sets, of which only one is used in the fit. Con- 
vergence is achieved when the quality of the fit to data which are not used for minimiztion 
stops improving. 

An important feature of our approach is that it is possible to check quantitatively the 
statistical features of results using suitable estimators. For example, one can check that the 
results do not depend on choices made during the fitting procedure, such as the choice of 
architecture of neural networks, which is analogous to the choice of parton parametrization 
in conventional fits. Namely, we repeat the fit with a different choice, and compute the 
distance 



d[q] 
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where q[ , q\ 2 ^ are the predictions for the i-th data point in the two fits, and c| , <t| 
the predictions for the corresponding statistical uncertainties, and the average is performed 
over all data. The results of the first and second fit are the same if d[q] = 1 on average. 
This also checks that the statistical uncertainties are correctly estimated. One can similarly 
check stability of the uncertainty estimate. In Ref. [5] this comparison has been performed 
succesfully. 

In Fig. [1] we compare our results for the NS structure function to other published 
determinations. These results are available through the webpage of the NNPDF Collabo- 



ration: http://sophia.ecm.ub.es/nnpdf. The large uncertainty that we find is a genuine 



feature of the determination of the nonsinglet quark distribution from the data included 
in our fit, and, especially at small x, it appears to reflect the current knowledge of the 
nonsinglet quark distribution. Indeed, for x < 0.05 the only data which constrain the (?ns 
combination in global fits are the data used in the determination of Ref. [S]. Hence, our 
results suggest that standard fits might be underestimating PDF uncertainties. 

In recent work on PDF uncertainties [H[S] it has been suggested that, mostly because of 
inconsistencies between data, the variation of the total \ 2 which corresponds to a one-sigma 
variation of the underlying PDFs is of order of A\ 2 ~ 50 for the global fits presented in 
those references instead of A% 2 = 1 of a statistically consistent fit [3] . In our approach, this 
quantity can be computed. We get A\ 2 ~ 1.7 (preliminary). This implies that the NMC 
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13 GeV 2 <Q 2 < 17GeV z 




Figure 1: 

The nonsinglet structure function as determined by the NNPDF collaboration [9] 
from 229 NMC and 254 BCDMS data points, compared to data and various other 

determinations. 

and BCDMS data are mostly consistent, though some inconsistent data are present [H1IH]- 
An extensive discussion of the way the published [5] and forthcoming [TU] fits based on the 
neural network approach can be used for the determination of physical parameters (such as 
a s ) and statistical properties of the data (such as A% 2 ) will be presented in a forthcoming 
publication. 

In [9] the strong coupling a s (M|) was fixed, but we could also extract it from the fit. The 
results of a preliminary analysis, shown in Fig. [21 suggest that nonsinglet data determine 
a s (M§) with an uncertainty which is rather larger than that (Aa g (M§) ~ 0.002) obtained 
in comparable determinations (see e.g. Ref. [llj). This preliminary result is consistent with 
that obtained using the same data in Ref. [12], with a method which eliminates the need 
to choose a parton parametrization. This strengthens the conclusion that uncertainties in 
available PDF fits might be underestimated. 

3 Towards a full parton set 

The extension of the results described in Ref. [5] to a full global PDF fit has benefited from 
the increased manpower of the NNPDF Collaboration, and is at a rather advanced stage [TU] , 
In particular, the evolution formalism of Ref. [9] has been extended to the computation of 
a full set of neutral-current and charged-current structure functions and fully benchmarked. 
A first full neural parton fit is in is in preparation. It will at first be based on DIS data only, 
including all available F% and fixed target data and the full NC and CC HERA reduced 
cross sections. 
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Figure 2: 

The x 2 profile for a preliminary NNLO determination of a s (M§ ) from NS data. The 
number of data points included in the fit is iVdat = 483. 

References 

[1] W. Giele et al. The QCD/SM working group: Summary report, hep-ph/0204316 . 
[2] M. Dittmar et al. Parton distributions: Summary report, ihep-ph/051 1119] . 

[3] Sergey Alekhin. Parton distributions from deep-inelastic scattering data. Phys. Rev., D68:014002, 2003. 

[4] J. Pumplin et al. New generation of parton distributions with uncertainties from global QCD analysis. 
JHEP, 07:012, 2002. 

[5] A. D. Martin, R. G. Roberts, W. J. Stirling, and R. S. Thornc. Uncertainties of predictions from parton 
distributions. I: Experimental errors. Eur. Phys. J., 028:455-473, 2003. 

[6] Stefano Forte, Lluis Garrido, Jose I. Latorre, and Andrea Piccione. Neural network paramctrization of 
deep-inelastic structure functions. JHEP, 05:062, 2002. 

[7] W. T. Giele, S. A. Keller and D. A. Kosower. Parton distribution function uncertainties. 
|hep-ph/0l04052| . 

[8] Luigi Del Dcbbio, Stefano Forte, Jose I. Latorre, Andrea Piccione, and Joan Rojo. Unbiased deter- 
mination of the proton structure function with faithful uncertainty estimation. JHEP, 03:080, 
2005. 

[9] Luigi Del Dcbbio, Stefano Forte, Jose I. Latorre, Andrea Piccione, and Juan Rojo. Neural network 
determination of parton distributions: The nonsinglet case. JHEP, 03:039, 2007. 

[10] Richard D. Ball, Luigi Del Debbio, Stefano Forte, Alberto Guffanti, Jose I. Latorre, Andrea Piccione, 
Juan Rojo, and Maria Ubiali. Neural network determination of parton distributions: The singlet case. 
In preparation, 2007. 

[11] J. Blumlein, H. Bottcher and A. Guffanti. Non-singlct QCD analysis of deep inelastic world data at 
0(alpha(s)**3). Nucl. Phys. B 774 (2007) 182 

[12] Stefano Forte, Jose I. Latorre, Lorenzo Magnea, and Andrea Piccione. Determination of ct B from scaling 
violations of truncated moments of structure functions. Nucl. Phys., B643:477-500, 2002. 



DIS2007 











PRELIMINARY RESULTS 




1 







